Evolution of a dense neutrino gas in matter and electromagnetic field 
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We describe the system of massive Weyl fields propagating in a background matter and interacting 
with an external electromagnetic field. The interaction with an electromagnetic field is due to the 
presence of anomalous magnetic moments. To canonically quantize this system first we develop the 
classical field theory treatment of Weyl spinors in frames of the Hamilton formalism which accounts 
for the external fields. Then, on the basis of the exact solution of the wave equation for a massive 
Weyl field in a background matter we obtain the effective Hamiltonian for the description of spin- 
flavor oscillations of Majorana neutrinos in matter and a magnetic field. Finally, we incorporate in 
our analysis the neutrino self-interaction which is essential when the neutrino density is sufficiently 
high. We also discuss the applicability of our results for the studies of collective effects in spin-flavor 
oscillations of supernova neutrinos in a dense matter and a strong magnetic field. 
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I. INTRODUCTION 

It is known that the neutrino interaction with exter- 
nal fields can significantly influence the evolution of su- 
pernova neutrinos For example, the combined ac- 
tion of a background matter and a magnetic field can 
cause the resonant transition like, v~ <-> vt (see, e.g., 
Ref. 0), where the indexes ± denote different helicity 
states. Hence the active neutrinos of the flavor a can be 
converted into sterile neutrinos of another flavor /3. 

Besides external fields, other factors, like the neutrino 
self-interaction, can strongly influence neutrino oscilla- 
tions. It happens, e.g., in a supernova explosion, when 
the typical neutrino luminosity can be ~ 10 51 erg/s 0. 
Accounting for the average supernova neutrino energy 
E v ~ lOMeV, we get that at the distance about several 
tens of kilometers from the protoneutron star surface the 
neutrino number density still can be high enough for in- 
teractions between neutrinos to be as important as the 
neutrino interaction with external fields. This neutrino 
self-interaction leads to the collective effects in neutrino 
oscillations. 

For the first time the neutrino self-interaction was con- 
sidered in Ref. [|| . Since then a lot of works on this sub- 
ject has been published (see, e.g., the recent review [5[ 
and references therein). Note that in the majority of the 
studies of collective effects in neutrino oscillations only 
the combination of the interaction with a background 
matter and the neutrino self-interaction was considered 
(see, e.g., Ref. Q). In the present work we shall general- 
ize the previous approaches for the description of collec- 
tive neutrino oscillations to include the interaction with 



an external electromagnetic field since, as we mentioned 
above, the influence of a strong magnetic field on the 
neutrino system evolution can be also important. 

Neutrinos can interact with an external electromag- 
netic field due to the presence of anomalous magnetic mo- 
ments. Note that the structure of the magnetic moments 
is completely different for Dirac and Majorana neutrinos 
(see, e.g., Ref. \J\). Despite the fact that nowadays there 
is no universally recognized confirmation of the nature 
of neutrinos in the present work we shall suppose 
that neutrinos are Majorana particles. Note that in var- 
ious scenarios for the generation of elementary particles 
masses, it is predicted that neutrinos should acquire Ma- 
jorana masses 

In the present work we shall resolve several important 
problems for the physics of Majorana neutrinos interact- 
ing with external fields. First, in Sec. lU using the re- 
sults of our recent paper [l(| we propose the classical field 
theory treatment of massive Weyl fields propagating in 
background matter and interacting with an external elec- 
tromagnetic field. Then, in Sec. Mil on the basis of the 
exact solution of the wave equation for Weyl fields in a 
background matter we canonically quantize these fields. 
In Sec. IIV1 in frames of our method we re-derive the 
effective Hamiltonian for the description of spin-flavor 
oscillations of Majorana neutrinos in matter and a mag- 
netic field. Finally, in Sec. [V] we apply the developed 
formalism to get the contribution of the neutrino self- 
interaction to the effective Hamiltonian. In Sec. IVI1 we 
summarize our results. 



II. CLASSICAL FIELD THEORY 
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In this section we develop the classical field theory de- 
scription of the massive neutrinos eigenstates, which are 



2 



supposed to be Majorana particles, in matter and elec- 
tromagnetic field. For this purpose we derive the Hamil- 
tonian for the system of two-component Weyl spinors and 
show that the classical canonical equations are equivalent 
to the wave equation for a Majorana neutrino in external 
fields. 

The wave equation for the neutrino mass eigenstates 
tp a , propagating in background matter and interacting 
with an external electromagnetic field, is known to have 
the form, 

[i-fB^-maWa- ^<W^6+ 5^7^6 = 0, (2.1) 

where m a are the masses of the particles, 7 M , 7 s , and 
&nv = (i/2)(7/i7^ — lulu) are the Dirac matrices. Note 
that we will formulate the dynamics of the system (|2.1[) 
in the mass eigenstates basis rather than in the flavor 
basis, as it is usually done when neutrino oscillations are 
considered, since only in the mass eigenstates basis one 
can distinguish between Dirac and Majorana masses fill ) . 

The interaction with matter is characterized by the ex- 
ternal fields (g^ b ), which, in principle, are nondiagonal in 
the neutrino mass eigenstates basis. In general case the 
matrix (g^ b ) is hermitian. However we shall discuss the 
situation when the CP invariance is conserved. Despite a 
current attempt to detect CP violating terms in the neu- 
trino sector [12j . no definite results have been obtained 
yet. In this case the vacuum mixing matrix is orthogonal 
and the matrix (g^ b ) is symmetric. The zero component 
of this matrix, (g® b ), contains the effective potentials of 
the neutrino interaction with non-moving and unpolar- 
ized matter, whereas the vector components, (g a b), are 
the linear combinations of the averaged matter velocity 
and the polarization. The explicit form of these matrices 
and the details of the statistical averaging can be found 
in Ref. 0. 

Note that the vector term in the neutrino matter in- 
teraction ~ g ablush is omitted in Eq. (|2.1j) since it is 
washed out for Majorana neutrinos. The contribution 
of the axial- vector interaction with matter to the wave 
equation (|2.1[) ~ g^bl^l 5 ^ is twice the analogous con- 
tribution for Dirac particles since both neutrinos and an- 
tineutrinos eq ually interact with a background matter 
(see, e.g., Ref. [i"4f). 

Neutrinos can interact with an external electromag- 
netic field F^ u — (E, B) owing to the presence of the 



anomalous magnetic moments (//„;,). It is know (see. e.g., 
Ref. [Hi p that the matrix (fiab) should be hermitian and 
pure imaginary, i.e. fi ab = -fi ba and \i* ab = -fi a b- We 
shall discuss the situation when no admixture of sterile 
neutrinos is in the mass eigenstates ij) a . In this case the 
electric dipole moments are equal to zero [15[ . 

We define the interaction with external fields in the 
mass eigenstates basis. However the interaction with a 
background matter is usually given for flavor neutrinos 
(see, e.g., Ref. 0). For the detailed discussion of the 
explicit forms of the matrices (g^ b ) and (/i Q fc) in the fla- 
vor eigenstates basis the reader is referred to the recent 
review [16[ . 

Since the neutrino mass eigenstates ip a are supposed 
to be Majorana particles they should obey the Majorana 
condition in the form, ip° = ij 2 ipa = ^cipa, where x c is 
a phase factor which we shall take equal to one. If we 
express the four-component Majorana spinors in terms 
of two-component Weyl fields, rj a and £ , as 

*<> - ('7") , or # . (_£. ) , (2,) 

which satisfy the Majorana condition, we can rewrite 
Eq. (|2.1I) in the two equivalent forms, 

?) a -(crV)ry a + m a a 2 i]* l - ^abc(B - iE)a 2 r] b 

+ K9ab + <rgab)Vb = 0, (2.3) 

or 

£a+(<rV)£ - m a a 2 C + Ma6<r(B + iE)a 2 4* 

" Kdab - °"*g«6)6 = 0. (2.4) 

Here er = (a\, a 2 , 03) are the Pauli matrices. In the fol- 
lowing we shall postulate these equations. Note that the 
analog of Eq. (12.3[) was previously derived in Ref. (l7j . 
We shall choose the spinors r\ a as the basic ones since it 
was experimentally confirmed that active neutrinos cor- 
respond to left-handed particles. Note that the rigorous 
proof of the equivalence of Majorana and Weyl fields was 
given in Ref. |18| . 

In Ref. [l(| we demonstrated that the classical dynam- 
ics of a massive Weyl field in vacuum should be described 
only in frames of the Hamilton formalism. Generalizing 
the results of Ref. [l(| to include the interaction with a 
background matter and an electromagnetic field we arrive 
to the following Hamiltonian: 



+ {^[tJ <r(B - iE)<7 2% * + r?J(720-(B + iEK] - i[nj{g ab + cr gab )ry b - (^) T ( 5 ° b + crg Qb )<] }] , (2.5) 

ab 

where 7r a are the canonical momenta conjugate to the "coordinates" r\ a . Using the aforementioned properties of the 
matrices (fJ, a b) and (g^ b ) we find that the functional (I2.5[) is real as it should be for a classical Hamiltonian. 
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Applying the field theory version of the canonical equa- 
tions to the Hamiltonian H, 

S H , 

i]a =t — = (<tVJ?7q - m a <j 2 r) a 

+ /i a &<r(B - iE)<j 2 Vb ~ Kdab + <TEab)rib, (2.6) 
SH 

TTa = - T = (<T V)7T a + m a <T 2 n a 

or)a 

- fi ab <J 2 (T{B + iE)< + i{g° ab + cr*g Qb )7r b , (2.7) 



one can see that in Eq. (|2.6j) we reproduce Eq. f|2 . 3[) for 
Weyl particles, which correspond to left-handed neutri- 
nos, interacting with matter and electromagnetic field. If 
we introduce the new variable £ a = io27r a , we can show 
that Eq. (|2.7|) is equivalent to Eq. (|2.4p for right-handed 
neutrinos. 

Previous quantum field theory based studies of Majo- 
rana neutrinos in an electromagnetic field and in a back- 
ground matter [2(| [Hj] involved the Lagrange formalism. 
The mass term in a Lagrangian for the Weyl field rj a has 
the form (see, e.g., Ref. (l8j). 



(2- 



It is however clear that C m vanishes if r\ a is a first- 
quantized field having commuting c- number components. 
To resolve this problem in Refs. |20L l2l| it was supposed 
that two-component Weyl fields are represented via the 
anti-commuting operators. Thus the treatment of Ma- 
jorana particles in those works can be considered just 
as the re-expression of already quantized fields in terms 
of the new variables rather than the generic canonical 
quantization. Moreover in Ref. [20( it was claimed that 
"there is no 'first-quantized' description of a massive two- 
component field in terms of c-number wave functions" . 

On the contrary, in the present work we have demon- 
strated that classical (first-quantized) massive Weyl fields 
in presence of an external electromagnetic field and a 
background matter can be perfectly described within the 
Hamilton formalism, cf. Eqs. (|2.5I) - (12.7[) . This fact just 
means that the Lagrange formalism is not a suitable tool 
for the studies of Majorana particles. Note that a more 
detailed description of the Weyl fields dynamics in vac- 
uum in frames of the canonical approach as well as the 
discussion of the applicability of the Lagrange and the 
Hamilton formalisms is presented in out recent work [Io| . 



III. QUANTIZATION 

In this section we canonically quantize a Weyl field 
propagating in a background matter. On the basis of an 
exact solution of the wave equation for a massive Weyl 
field in matter we express the energy and the momen- 
tum of the field as a sum of contributions of independent 
quantum oscillators. From the requirement of the posi- 
tive definiteness of the energy it turns out that the oper- 



ators in the decomposition of the wave functions should 
obey the Fermi-Dirac statistics. 

In the following we shall suppose that the background 
matter in average is at rest and unpolarized, i.e. g ab = 0. 
This approximation is valid in almost all realistic cases. 
Indeed the matter motion is relevant for the neutrino 
dynamics if the speed of medium is comparable with the 
speed of light. This situation may be implemented, e.g., 
if a beam of neutrinos propagates inside a relativistic jet 
from a quasar. Since the results of the present work are 
likely to be applied for the studies of supernova neutrinos, 
the matter motion seems to be irrelevant for us. 

When neutrinos interact with a non-degenerate 
plasma, we may neglect the matter polarization if 
/i/|B| <C T (see Ref. [19] ) , where [if ~ is the magnetic 
moment of a background fermion, /ie is the Bohr mag- 
neton, and T is the plasma temperature. In the present 
work we shall study the influence of both external fields 
(see Sec. lIV[) and the neutrino self-interaction (see Sec. lV)) 
on neutrino oscillations. In Ref. (22[ it was was found that 
the collective effects in neutrino oscillations reveal them- 
selves most intensively at the distance r ~ 100 km from 
a protoneutron star. The magnetic field at this distance 
can be B = B (R/rf ~ 10 10 G, where B ~ 10 13 G is 
the typical magnetic field on the protoneutron star sur- 
face and R ~ 10 km is the stellar radius. Supposing 
that T ~ 1 MeV [23| , we get that the matter polariza- 
tion becomes unimportant for stars possessing moderate 
magnetic fields. 

Let us decompose the Hamiltonian (I2.5[) into two terms 
H = Hq + Hi nt . The former term, Hq, contains the 
vacuum Hamiltonian as well as the matter interaction 
term diagonal in the neutrino types, 



H = f d 3 r^{7r Q T [(<TV)-i.g Q Q K 

-K) T [(-v)-i 9 2>: 

+ m a [{V*a) T V2Ka + « ) T CT 2 ?7a] }■ 



(3.1) 



The latter term in this decomposition, 



H, 



f d 3 r^{/i ab [7rJ CT (B-iE)(7 2 ry 6 * 

J a^b 

+ j7j(7 20 -(B + iEK] 
-ifab [njr]b - (?7a) T <] }> 



(3.2) 



has the nondiagonal matter interaction and the interac- 
tion with an electromagnetic field which is nondiagonal 
by definition. 

Analogously to Eqs. (|2 .6[) and (|2.7p we define the re- 
duced Hamilton equations which contain only the Hamil- 
tonian H : r]^ = 5Ho/5iVa and 7ri 0) = -6H /6rj^K 
Using the results of Refs. [13, [13 we can find the solu- 
tions of these equations in the form, 
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^ 0) M) 



£i 0) M) =i 



d 3 p 

(2tt) 3 / 2 



a n W-e 



Et + |p| - 9. 



E a + |p| + g° aa 



\E~t 



d 3 p 

(2tt) 3 / 2 



l+ -iE+t , TO 



-ipr 



(&-)*^ + e iB «< 



+|p|+s2a a 

^ (b + )*w e iEit 



(3.3) 



where we introduce the new variable ^i ^ = io^i ' 1 , w± are the helicity amplitudes defined in Ref. [24j . and 



4° = V™2 + (IpI-C<7°J 2 



(3.4) 



is the energy of a Weyl field [T3, H3, C — ±1 is the particle helicity. To derive Eqs. p.3[) and Q3.4p we suppose that 
the external field g® a is spatially constant. 

Using Eq. (|3.3p we can express the Hamiltonian (J3TTJ) as 



Ho=j d" P J2{E- [K(p)]*6-(p) + [&"(p)]*a-(p) 



+ 



Pfr o ^ 2 i a ° (p) (p)r + b « (p) [ fl « (p)] * ) 



(£« + | P |+S°J 



a^(p)^(p)]*+6^p)[a+(p)]' 



(£+ + |p| - 5aa 



f— — {[6+(p)]*a+(p) + [a a +(p)]*6^(p)} 



la { 



#a + IPl + ffga 



[e- 2i ^* {o-(p)6-(-p) + b-(-p)a-(p)} 



+ e 2i ^*{[a-(p)]*[&-(-p)]* + [6-(-pF[ a -(p)]*}] 



+ 



E+ 



E+ + \p\-g 



- [e- 2iE « 1 {a+(p)6+(-p) + 6+(-p)a+(p)} 



2i£ « +t {[a+(p)]*[6+(-p)]* + [6+(-p)]*[a+(p)r}]}}, 



(3.5) 



in terms of the creation, [ffl„(p)]* and [6^ (p)]*, as well 
as the annihilation, a„(p) and 6„ (p), operators. Note 
that the operators a„ (p) and 6„ (p) are independent up 
to now. 

Let us establish the following relation: 
a±(p)(E± + | P | TflSa) = 46„ (p)(|p| T5q Q ). (3-6) 

and the analogous expression for conjugate operators, as 
well as suggest that the operators a„ (p) obey the anti- 
commutation properties, 

{a±(k);[a±(p)n + =^ 3 (k-p), (3.7) 

with all the rest of the anticommutators being equal to 
zero. In this case the time dependent terms in Eq. (13. 5[) 
are washed out. The remaining terms can be represented 



H =[ d 3 p J2iE- (a- )* a- +E+(a+y at] 

+ divergent terms, (3-8) 

which shows that the total energy of a massive Weyl field 
is a sum of energies corresponding to elementary oscilla- 
tors of positive and negative helicities. 

Using the results of Ref. we can also quantize the 
total momentum of a Weyl field defined as 



J o=/d 3 r£[(^ 0) *) Tv ^ 

J a L 

-(-i 0) ) T v^ 



(3.9) 
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With help of Eqs. (j3^)l . and (J3J]) we rewrite 

Eq. p. 91) in the following form: 

P = / d 3 p^p[K-)*a- + (a+)*a+] 

+ divergent terms, (3.10) 

which has the similar structure as Eq. (|3.8|) . 

Note that the divergent terms in Eqs. (|3.8[) and (13. 10)) 
contain the factor <5 3 (p = 0) — > oo, which can be formally 
removed by the normal ordering of the operators and 
(a„)*. It is also interesting to mention that a massive 
Weyl field in vacuum can be quantized in the two inde- 
pendent ways (see Ref. [l(|) because of the degeneracy 
of the neutrino energy levels: E~ = E+ = \/m^ + |p| 2 . 
On the contrary, in matter only one of the possibilities 
for the quantization gives the correct result for the total 
energy (|3.8p since the energy levels are no longer degen- 
erate, cf. Eq. (EH)) . 



Note that Eqs. (|4.1[) - (|4.3I) are valid for arbitrary neutrino 
masses, momentum, and the diagonal neutrino interac- 
tion with matter. 

Now we define the neutrino density matrix as 

<S 3 (p - k) PAB (p) = (a^(p)o A (k)>, (4.4) 

where A = (C,a) is a composite index and (...) is the 
statistical averaging over the neutrino ensemble. In prin- 
ciple, we could interchange the indexes A and B in the rhs 
of Eq. I|4.4p . Since we study the dynamics of Majorana 
neutrinos, such a transposition would mean the consider- 
ation of neutrinos as antiparticles rather than as particles 



IV. NONDIAGONAL INTERACTION WITH 
MATTER AND ELECTROMAGNETIC FIELD 

In this section we apply the approach for the quantiza- 
tion of a massive Weyl in a background matter developed 
in Sec. IIIII for the treatment of the nondiagonal Hamil- 
tonian H- lnt given in Eq. (|3.2[) . On the basis of the ob- 
tained results and using the density matrix formalism (25j 
we derive the effective Hamiltonian for the description of 
neutrino spin-flavor oscillations in matter and electro- 
magnetic field. Then we demonstrate that for ultrarel- 
ativistic particles our effective Hamiltonian is consistent 
with the analogous expression obtained in frames of the 
standard quantum mechanical approach. 

To quantize the Hamiltonian H lni we shall use the for- 
ward scattering approximation. It means that one has to 
account for only the terms conserving the number of par- 
ticles p(i . Using Eqs. ([33]), dH]), and ([3J]) we rewrite 
Eq. (pOj) in the form, 



(4.1) 



(4.2) 



(4.3) 

I 

as we do here. Anyway both definitions will give equiv- 
alent results. Note that a density matrix, having both 
neutrino type and helicity indexes, which is analogous to 
Eq. (|4.4[) . was studied in Ref. [13] where the interaction 
between Dirac neutrinos, mediated by a scalar boson, was 
discussed. 

Applying the quantum Liouville equation for the de- 
scription of the density matrix evolution, 

ip=\p,H- mt ], (4.5) 

we can rewrite it as ip = [H, p] using the effective quan- 



ffmt = J d 3 P ^[Af Q -(a-)* a - e i5 - t + M+(a+ya+^ t + F^a~ )* a+ + F+(atT^e^% 



where <5± = £± - E±, <j± = £± - Ef, 



and 



M ± = _ 5ob ( E£ + \p\T9°aa , Ef + |p| T 9° bb 
ab + 4 ^ | P |T.9 Q Q | P |T.9 b ° b 



m a m b 



(K ± + | P |T5 a ° a )(^ ± + |p|T3 fc ° b ) 



± 



Mafe j 1 

4 \|p|T52a 

1 



m a wl(cT[B =F iE])< + m b ™ T (<x[B ± iE])w* S +l , Pl , Tg ° Q a 

Et + IpI T 9° bb 



|P|T5 b u b 



m b w;^(cr[B ± iE])w+ + TO a w T (cr[B T iE])i 



^ + |p|T.g b ° b 

E± + \p\ T g° aa 



E± + \p\T9° aa , ET + \ p \±g° b 



\p\±9° bb 



i a m b wZ(cr[B =F iE])u>± 



{Et + |p| T9° aa )(\Pl ±9° bb ) (E^ + |p|±«&)(|p| T9° ao 
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turn mechanical Hamiltonian, 



H 



(M-e^) (F-e 




(4.6) 



Again we stress that yet no expansion over the parame- 
ter m a /|p|, which is small for ultrarelativistic neutrinos, 
is made. Thus Eqs. (|4.5|) and (|4. 6[) are valid for the de- 
scription of neutrinos with arbitrary initial momentum. 

To study the evolution of our system we use Eq. (|3.3p , 
where the wave functions already contain time dependent 
exponential factors. The energies in Eq. ()3.3|) correspond 
to the total diagonal Hamiltonian (|3.1[) . cf. Eqs. Q3.4p 
and (|3.8I) , which contains both the mass term and the di- 
agonal interaction with a background matter rather than 
only a kinetic term as in Ref. [25| . Thus our treatment 
is analogous to the Dirac picture of the quantum theory. 
That is why in Eq. (|4.5p it is sufficient to commute the 
density matrix only with H mt rather than with the total 
Hamiltonian H = Hq + Hi nt . 

To derive Eqs. (|4.1j) - (|4.3j) we suppose that the external 
fields g° b , a ^ b, and F^ are spatially constant. In fact, 
for this supposition to be valid, the characteristic scale 
of the external field variation L cxt should be much bigger 



than the typical width of the neutrino wave packet h/ E v . 
One can see that in almost all realistic situations external 
fields can be regarded as constant for the derivation of an 
effective Hamiltonian. Nevertheless, when the obtained 
effective Hamiltonian is used to describe the dynamics of 
the system, e.g., with help of Eq. (|4.5|) . we should account 
for the variation of external fields. 

To demonstrate the consistency of the obtained results 
with the achievements of the standard quantum mechan- 
ical description of neutrino spin-flavor oscillations (see, 
e.g., Ref. @]) we discuss the simplest case of the two neu- 
trino eigenstates, a — 1,2, and consider the situation of 
ultrarelativistic particles, |k| ^> max(m , <7° a ), where k 
is the initial momentum of neutrinos. In this case we 
should decompose the energy levels (|3.4p as 



E; 



2|ki 



T9° aa 



(4.7) 



Then, supposing that E = 0, since it is difficult to create 
a large scale electric field, we get that M^ b ss T9ab and 
F^ b ss — /j, ab | B | sin^kB, where $ k B is the angle between 
the vectors k and B. 

To eliminate the time dependent factors in the effec- 
tive Hamiltonian (I4.6[) we make the transformation of the 
density matrix 17], 



Pqm = UpU\ U = diagle- 1 ^^)*, S*-&*, e -i(»-*?i)*, «.«*+«&)*}, (4.8) 

where $ = <5m 2 /4|k| is the phase of vacuum oscillations and 8m 2 — m\ — m\ is the mass squared difference. 

The evolution of the transformed density matrix can be represented as ip qm = [H qm , p qm ], where the new effective 
Hamiltonian has the form, 

/ $ + 9i2 -^i 2 |B|sintf k B\ 

% ea =UMft + JlM=\ 9 k + -^i|B|sintf kB (49) 

q -^ 12 |B|sini?kB -.9i2 

\-M 2 i|B|sin^ kB -g° 21 -$ - g% 2 J 



Recalling the properties of the magnetic moments ma- 
trix for Majorana neutrinos, pn = ifJ- and (121 = — 
with pi being a real number, we can see that Eq. (|4.9p 
reproduces the well known quantum mechanical Hamil- 
tonian for spin-flavor oscillations of Majorana neutrinos 
in matter and a magnetic field 0- 

Note that previously calculated transition probabilities 
of neutrino oscillations in a magnetic field [2fJ| and in a 
background matter [2l| can be re-derived using the effec- 
tive Hamiltonian (|4.9p which was obtained in frames of 
the approach involving canonically quantized Weyl fields. 
However, as we mentioned in Sec. [Hi methodologically 
our method for the description of massive Majorana neu- 
trinos in external fields is more logical since it is based 
on the first principles of the quantum field theory. 



V. SELF-INTERACTION 

In this section we generalize the results of Sees. HB 
IIVI to include the neutrino self-interaction. First we re- 
formulate the previously proposed Hamiltonian for the 
self-interaction in terms of the two-component Weyl fields 
and then we quantize it. Finally, using the density matrix 
formalism we derive the corresponding contribution to 
the quantum mechanical effective Hamiltonian and com- 
pare it with the previously obtained results. 

The Hamiltonian describing the neutrino self- 
interaction, mediated by a neutral Z-boson, was derived 
in Refs. @, [25| and has the form, 

H s = J d 3 rJ2 GabG^^b-^^d, (5.1) 

abed 
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where 7^ = 7^(1 — 7 5 )/2 and G a b are the coeffi- 
cients which depend on the neutrino interactions chan- 
nel. The explicit form of these coefficients can be found in 



Ref . 28j . Let us choose the wave function tf> a in Eq. (|5.ip 



as ijj^ in Eq. (|2 . 2|) . Therefore we express the Hamilto- 
nian (|5.1|) in terms of the two-component Weyl spinors 
as 



H s = J d 3 r GabGcdTjla^b ■ r)t^r] d , (5.2) 

abed 

where ct m = (1,<t). It should be noted that we presented the heuristic derivation of Eq. (|5.2[) from Eq. (|5.ip . In 
principle we could just postulate Eq. (|5.2I) . 

Using Eqs. (13.31) and (I3.6[) we cast the self-interaction Hamiltonian (|5.2[) into the form, 

H S = -^3 J d 3 pd 3 p'd 3 qd 3 q'<5(p + q - p' - q') £ G ab G cd [7 h ( q )*} T a^') ■ Mp)*} T ^Vd(p'), (5.3) 

^ ^' abed 

where 

a(p) = £ [°i C) (p)«i C) (p)e- ,B - 0t + [4 C) (-P)f^ , (-P)^ ,t ], (5.4) 



Va 

C=±i 



is the Fourier transform of the wave function 77 a and 

\ (P) = W-(P), (P) = w_(p), u+(p) = - -=+ — j-j— w+(p), ?j-(p) = — j— — — Q w+(p), (5.5) 

+ |P| Saa + |P| + 9, 



aa 



are the basis spinors rewritten in a formalized manner, cf. Eq. (|3.3p . 

Applying Eq. (|4.5[) to account for the contribution of the self-interaction to the dynamics of the system and again 
working in the forward scattering limit p6l | , we can represent the evolution of the density matrix, ip = [Hs , p] , using 
the effective Hamiltonian Hs ■ After a bit lengthy but straightforward calculations we get the following expression for 

n s -. 

/i3 
^{M^(k,k)tr{[M^(p,p)-iV^(p,p)]p(p)} + ^(k,k)tr{[^(p,p)-M' 1 (p,p)]p(p)} (5.6) 

- [M M (k,p)-JV M (k 1 p)]p(p)[M"(p,k)-JV(p J k)] + (if l ,(k J p) - [^(k,p)] T V T ( P )(L^(p,k) - [L^( P ,k)] T )}, 

where 

M^ B (p,k) =G ab e iE ^ t (u A (p)\a»\u B (k)}e- iEB W t , i^ B (p,k) = Gabe^^iuAip^lvsik))^^^, 
N» B (p,k) =G ba e- iEB ^ t {v B (k)\a^\u A (p))e iEA ^ t , L AB (p,k) = G afce - iB ^< : (7u(p)|aHMk)> e - ifi «( k )' . (5.7) 

In Eq. (|5.6p the transposition means the interchange of both discrete and continuous indexes, i.e. [L AB (p, k)] T = 
G 6a e- iiSf3(k)t (7;B(k)| ( 7^|M A (p))e- i ^(p)* etc. 

Note that Eqs. (15 .6[) and (|5.7I) are valid for arbitrary neutrino masses, initial momenta, and the diagonal neutrino 
interaction with matter. However the analysis of these expressions is quite cumbersome. That is why again we 
discuss the situation of the two neutrino generations, a = 1,2, and suppose that neutrinos are ultrarelativistic 
particles, |k| 3> max(777 a , g^ a ). Then, to eliminate the time dependence in Eq. (|5.7p we make the additional matrix 
transformation of the effective Hamiltonian, TL qm — UHs^ , where Hs — Hs[MpM' f ], since %s is the function of p, 
and the matrix U is defined in Eq. (|4.8p . 

Finally we can represent the contribution of the neutrino self-interaction to the quantum mechanical effective 
Hamiltonian as, 

H qm (k) = diag(^__,H++), (5.8) 

where 

U — =2 / " C0S ^ P ){Gtr[Gp__(p) - G T P++ (p)] + G[ P __(p) - pT+(p)]G}, 
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n ++ =2 J - cos^ k p){G T tr[G T p ++ (p) - Gp__(p)] + G T [ P++ (p) - p T _(p)]G T }, (5.9) 

and TL± T = 0. Here is the angle between the vectors k and p and we use the helicity components of the density 
matrix, 



To derive Eqs. (15.81) and (I5.9|) we use the identity for 
two-component c-number spinors, 77 J ct m 772 • r]l a>J " r l'i — 
—r^jfj]^ ■ r/la^rj^, which results from the Fierz trans- 
formation of four-component spinors, ^17^(1 — J 5 )ip2 • 
"037^(1 - 7 5 )"04 = ~ 7 5 )V>4 • -027 M (1 - 7 5 )"03- 

One can conclude from Eq. (|5.9[) that the self- 
interaction influences spin-flavor oscillations of neutri- 
nos. However, since the nondiagonal terms in Eq. (15. 8[) 
are equal to zero for ultrarelativistic particles, the self- 
interaction cannot directly induce transitions between 
different helicity states. 

It should be noted that in the majority of works where 
collective effects in neutrino oscillations were studied, the 
case of Dirac neutrinos was examined. Although one can 
expect that in the ultrarelativistic case the dynamics of 
Dirac and Majorana neutrinos should be similar, we can- 
not reach a complete coincidence because Dirac particles 
have twice more degrees of freedom, i.e. an additional 
density matrix for antineutrinos is required. Neverthe- 
less let us check the consistency of our findings with the 
previously obtained results. First we should chose a def- 
inite helicity. For example, we may put H and 

W++ = 0. Then, defining the "antineutrino" density 
matrix as p = cf. Eq. (I4.4[) . we re-derive the con- 

tribution of the neutrino self-interaction to the effective 
Hamiltonian obtained in Ref. [25| . 

At the end of this section we notice that the pre- 
sented derivation of Eqs. (|5.8p and (|5.9I) is not unique. 
In the general self-interaction Hamiltonian (15.11) we can 

set i\) a = if)£\ where ipP is defined in Eq. (|2~2l . Thus 
the Hamiltonian H$ can be expressed in terms of the 
canonical momenta, 

H s = J d 3 r J2 GabGcd^a^t • Trja^, (5.11) 

abed 

where use the relation between £ a and 7r a : £ a = i<727r a 
(see also Sec. HI]) . 

Then we should re-define the wave functions r\ a and 
£ a in Eq. (|3.3[) , introducing the additional multiplier 1 /2 
in each spinor, as well as the connection between oper- 
ators a^(p) and 6„ (p), which now reads, a^{p)(E^ + 
\p\T9aa) = bt(p)(\p\ T9°aa): cf. Eq. Note that 

these modifications do not affect the results of Sees. IIIII 
and IIV1 Performing the same calculations which led us 
to Eqs. (|5.8|) and (|5.9p . but using the modified Hamil- 
tonian (|5.11l) . we get that for ultrarelativistic particles 



I 

the contribution of the neutrino self-interaction to the 
quantum mechanical effective Hamiltonian coincides with 
Eqs. (0 and 



VI. CONCLUSION 

In summary me mention that in the present work we 
have constructed the consistent quantum theory of a sys- 
tem of massive Weyl fields propagating in a background 
matter and interacting between themselves and with an 
external electromagnetic field. We have obtained several 
important results. 

First, in Sec. HH the classical field theory description 
of a massive Weyl field in an arbitrarily moving and po- 
larized matter and an electromagnetic field has been pre- 
sented. Using the approach developed in Ref. [To|, where 
the evolution of a massive Weyl field in vacuum was stud- 
ied, we have derived the classical Hamiltonian (|2.5p for 
our system. Then applying the canonical Hamilton equa- 
tions (|2.6[) and (|2.7[) we have re-obtained the analog of 
the well known wave equation (I2.1[) for a Majorana neu- 
trino in matter and an electromagnetic field. This our 
result corrects the previous statement [2(| that massive 
Majorana particles are essentially quantum objects de- 
scribed only using the creation and annihilation opera- 
tors formalism. Moreover, now we have expanded the 
classical field theory approach, cf. Ref. [l(J, to include 
the interaction with matter and an electromagnetic field. 

Second, in Sec. lIIH we have canonically quantized mas- 
sive Weyl fields in a non-moving and unpolarized matter. 
We have used the plane wave solution (|3.3[) and p. 41) (see 
also Ref. [13) of the corresponding wave equation, where 
we supposed that the expansion coefficients are the op- 
erators. Then, requiring the positive definiteness of the 
total energy (|3.5[) , we have obtained that the operator ex- 
pansion coefficients should satisfy the anticommutation 
properties Q3.7p . It is interesting to mention that unlike 
the quantization of a Weyl field in vacuum, where two 
independent quantization schemes are possible, in mat- 
ter there is only one opportunity (|3.6I) . which gives the 
correct form for the total energy (|3.8I) and the total mo- 
mentum (|3.10p . 

Third, in Sec. IIV1 we have applied the elaborated 
quantization method for the nondiagonal interaction with 
matter and an electromagnetic field p. 21) . Within the for- 
ward scattering approximation we have derived the quan- 
tized interaction Hamiltonian (|4.ip - (|4.3l) . which is valid 
for arbitrary neutrino masses, an initial momentum, and 
the diagonal neutrino interaction with matter. Then, us- 
ing the density matrix formalism, developed in Ref. |25| j , 
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and in the approximation of ultrarelativistic particles we 
have re-derived the effective Hamiltonian (|4.9p . previ- 
ously obtained in frames of the standard quantum me- 
chanical approach [2], for the description of spin- flavor 
oscillations of Majorana neutrinos in matter and a mag- 
netic field. 

Finally, in Sec. |V] we have quantized the self- 
interaction (15. ip of Majorana neutrinos using the devel- 
oped formalism. Again in the forward scattering ap- 
proximation we have got the contribution to the effec- 
tive Hamiltonian (|5.6[) and (j5.7p which is valid for neu- 
trinos with arbitrary masses and an initial momentum. 
Then, for ultrarelativistic particles we have compared our 
results with the previously obtained effective Hamilto- 
nian [25[, which describes collective neutrino oscillations, 
and have found the consistency. 

Note that in all the previous works where collective ef- 
fects in neutrino flavor oscillations were studied the case 
of Dirac neutrinos in a background matter was consid- 
ered. Thus in the present work for the first time we have 
discussed the situation of Majorana neutrinos and gener- 
alized the consideration to include an external magnetic 
field. It should be noted that for supernova neutrinos 
both the interaction with a dense background matter, a 
strong magnetic field, and the neutrino self-interaction 
can be of equal importance. Therefore, the effective 
Hamiltonians (14.91) , (|5.8|l , and (|5.9p , derived in our work, 
may be used for the treatment of collective effects in spin- 



flavor oscillations of supernova neutrinos in matter and 
a magnetic field. 

The wave functions Q3.3P exactly take into account the 
diagonal neutrino interaction with background matter, 
g® a , whereas in Sec. |IV] the external fields, g° b , a ^ 6, 
and F^ v are treated perturbatively, with only linear term 
being kept, cf. Eqs. (|4.2p and (|4.3p . However in general 
case the potentials <?° and g° b can be of the same order 
of magnitude. Nevertheless, besides supernova neutri- 
nos discussed here, we may apply the obtained results 
to study nonperturbative effects in neutrino oscillations 
if the nondiagonal potential of matter interaction g® b is 
small or negligible. For instance, such a situation is im- 
plemented when neutrinos propagate in the inner crust 
of a neutron star, where n e . p <C n n , with n e ^ p ^ n being 
the number densities of electrons, protons, and neutrons 
respectively. Another example, when the nondiagonal 
matter interaction is unimportant, is v p <-> v T oscilla- 
tions channel. 
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